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High energy heavy-ion collisions are proposed as a novel means to obtain in- 
formation about the high density (HD) behaviour of nuclear symmetry energy. 
Within an isospin-dependent hadronic transport model using phenomenological 
equations of state (EOS) for dense neutron-rich matter, it is shown that the 
isospin asymmetry of the HD nuclear matter formed in high energy heavy-ion 
collisions is determined mainly by the HD behaviour of nuclear symmetry en- 
ergy. Experimental signatures in several sensitive probes, i.e., 7r~ to 7r + ratio, 
transverse collective flow and its excitation function as well as neutron-proton 
differential flow, are investigated. A precursor of the possible isospin separa- 
tion instability in dense neutron-rich matter is predicted to appear as the local 
minima in the excitation functions of the transverse flow parameter for both 
neutrons and protons above the pion production threshold. Because of its quali- 
tative nature unlike other quantitative observables, this precursor can be used as 
a unique signature of the isospin dependence of the nuclear EOS. Measurements 
of these observables will provide the first terrestrial data to constrain stringently 
the HD behaviour of nuclear symmetry energy and thus also the EOS of dense 
neutron-rich matter. Implications of our findings to neutron star studies are also 
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discussed. 
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I. INTRODUCTION 



The rapid advance in nuclear reactions using rare isotopes has opened up several new 
frontiers in nuclear sciences • In particular, the intermediate energy heavy rare isotopes 
currently available at the National Superconducting Cyclotron Laboratory (NSCL/MSU) 
and the more energetic ones to be available at the future Rare Isotope Accelerator (RIA) in 
the United States provide a unique opportunity to explore novel properties of dense neutron- 
rich matter that was not in reach in terrestrial laboratories before. This exploration will 
reveal crucial information about the equation of state (EOS) of neutron-rich matter. To 
understand the latter and its astrophysical implications, such as, the origin of elements, 
structure of rare isotopes and properties of neutron stars are presently among the most 
important goals of nuclear sciences. The EOS of neutron-rich matter of isospin asymmetry 
5 = (p n — p p )/(p n + Pp) can be written as 

e(p,5) = e(p,0) + E sym (p)S 2 (1) 

within the parabolic approximation (see e.g., 0), where e(p, 0) is the energy per nucleon in 
isospin symmetric nuclear matter. To study the density dependence of nuclear symmetry 
energy E sym (p) has been a longstanding goal of extensive research with various microscopic 
and/or phenomenological models over the last few decades, e.g., |TD|-p2|, for a recent review, 



see, e.g. ||. All models are able to give the symmetry energy at normal nuclear matter 
density E sym (po) of about 34±4 MeV in agreement with that extracted from the atomic mass 
data. However, the predicted results on the density dependence, especially at high densities, 
are extremely diverse and often contradictory. At subnormal densities, all models predict 
that the symmetry energy increase with density although the growing rate is rather model 
dependent. To extract the density dependence of symmetry energy at low densities, several 
approaches have recently been proposed. These include measuring the thickness of neutron- 
skins of rare isotopes [^] and heavy stable nuclei, such as 208 Pb p3-27j, isospin fractionation 



Of as well as nuclear collective flow [31-33] in intermediate energy heavy-ion collisions. 
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Above normal nuclear matter densities, even the trend of the density dependence of E sym (p) 
is still controversial. Theoretical results can be roughly classified into two groups, i.e., a 
group where the E sym (p) rises monotonously and one in which it falls with the increasing 
density above about twice the normal nuclear matter density. The predictions also depend 
on the nature and form of effective interactions used. For instance, within the same Hartree- 
Fock approach using about 25 Skyrme and Gony effective interactions that have been widely 
used successfully in studying saturation properties of symmetric nuclear matter and nuclear 
structures near the (3 stability valley, the calculated symmetry energies were found to fall 



approximately equally into the two groups [24 ,34 



The density dependence of nuclear symmetry energy, especially at high densities |35| , has 



many profound consequences for various studies in astrophysics [^438[. In particular, an 



increasing E sym (p) leads to a relatively more proton-rich neutron star whereas a decreasing 
one would make the neutron star a pure neutron matter at high densities. Consequently, the 
chemical composition and cooling mechanisms of protoneutron stars p9|^0| , critical densities 



for Kaon condensations in dense stellar matter [EILO, mass-radius correlations |4"3|,p0 



as 



well as the possibility of a mixed quark-hadron phase H5[ in the cores of neutron stars will all 



be rather different. The high density (HD) behaviour of nuclear symmetry energy E sym (p) 
is very important for understanding many interesting astrophysical phenomena, but it also 
subjects to the worst uncertainty among all properties of dense nuclear matter fl35"f . The 
fundamental cause of the extremely uncertain HD behaviour of E sym (p) is the complete lack 
of terrestrial laboratory data to constrain directly the model predictions. We explore in this 
work the possibility of using high energy heavy-ion collisions to probe the HD behaviour 
of E sym (p). A brief report of this work can be found in ref. j46|. The E sym (p) is found 



to affect significantly several aspects of high energy heavy-ion collisions. In particular, the 
neutron/proton ratio of HD nuclear matter formed in these collisions is determined mainly 
by the HD behaviour of E sym (p). Several promising signatures of the HD behaviour of 
E sym (p) are also investigated. Besides the quantitative observables, we also study one special 
qualitative signature of the HD behaviour of the symmetry energy, i.e., the local minimum 



in the excitation function of nuclear transverse flow. The paper is organized as follows. In 
Section 2, we discuss the equation of state (EOS) for neutron-rich matter, in particular, the 
HD behaviour of nuclear symmetry and its effects on the neutron/proton ratio in neutron 
stars. In section 3, we study the neutron/proton ratio and its dependence on the symmetry 
energy in HD hadronic matter formed in high energy heavy-ion collisions within an isospin- 
dependent hadronic transport model. In section 4, we explore several experimental probes 
of the HD behaviour of nuclear symmetry energy. Finally, we summarize in section 5. 

II. EQUATION OF STATE OF DENSE NEUTRON-RICH MATTER 

In this section properties of the equation of state of dense neutron-rich matter are dis- 
cussed. In particular, we discuss the HD behaviour of nuclear symmetry energy and poten- 
tials as well as their influences on neutron stars. Based on these properties, we present our 
expectations about their influences on heavy-ion collisions involving neutron-rich nuclei. 



A. High density behaviour of nuclear symmetry energy 

In the parabolic approximation of the EOS for neutron-rich matter in Eq. [j], the sym- 
metry energy is [|L7 



E sym {p) = e(p, 1) - e(p, 0) = -E kin (p, 0) + V 2 (p), (2) 
where Eki n {p,0) is the kinetic energy per nucleon in symmetric nuclear matter 

and V^p) is the deviation of the interaction energy of pure neutron matter from that of 
symmetric nuclear matter. The E sym (p) becomes negative if the condition [|35| 

V 2 (p)<~E kin (p,0) (4) 
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is reached at high densities. A negative symmetry energy at high densities implies that the 
pure neutron matter becomes the most stable state leading to the onset of the isospin separa- 
tion instability (ISI). Consequently, pure neutron domains or neutron bubbles surrounding 



isolated protons may be formed in the cores of neutron stars f55fl . Energetic nuclear re- 
actions with rare isotopes provide a unique opportunity to pin down the E sym (p) at high 
densities and study the possible existence of ISI and its consequences. We use the following 
two representatives of the symmetry energy as predicted by many body theories 

KyM = E S ym(po)u (5) 

and 

EsymiP) = E sym(po)u • /_ , (6) 

where u = p/p and u c = p c J p is the reduced critical density at which the E b sym (p) crosses 
zero and becomes negative at higher densities. The predicted value of u c ranges from about 
2.7 (Hartree-Fock with the Skyrme interaction Sp j33J) to 9 (variational many-body approach 
with the UV14+UVII interaction flr7|.). For our numerical calculations in the following we 
use u c = 3. Our conclusions in this work are qualitatively independent of this value. Thus 
future comparisons with experiments may rule out the existence of ISI below this density, 
but they can not exclude the occurrence of ISI at higher densities. By using higher values for 
u c , effects found here will be simply shifted to higher beam energies. Therefore, excitation 
function studies in both theories and experiments will be necessary to determine the exact 
value of u c if it does exist in nature. The above two forms of E sym (p) are shown in the upper 
window of Fig. 1. By design, they both have the same value of E sym (p ) = 30 MeV at the 
normal nuclear matter density p an d are very close to each other at lower densities. At 
high densities they have completely different trends reflecting the diverging predictions of 
nuclear many-body theories. 

It is necessary to mention that in a recent publication P7 ], it was found that the symmetry 



enery itself is isospin-dependent at all densities. This amounts to absorbing the higher order 
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terms in 5 into the E sym in the expansion of the EOS of neutron- rich matter in Eq.(l). 
Because of the charge symmetry of nuclear interactions used the odd terms in 5 should 
vanish. The lowest order correction to the parabolic approximation in Eq. (1) is the <5 4 
term. The magnitude of the latter has been found extremely small in all microscopic many- 
body calculations (see, in particular those in refs. | ^T] , [I^[ , [I7] , |IB| , ^ , ^5| ) , except that in P7 



The origin of the unusual prediction in the latter deserves further investigations. In this 
work we stick to the parabolic approximation of the EOS of neutron-rich matter. 



B. Symmetry energy and proton fraction in neutron stars at (3 equilibrium 

The above two forms of the symmetry energy lead to significantly different predictions 
on several properties of neutron stars. To a good approximation suitable for this study, 
a neutron star without neutrino trappings can be considered as a npe-matter consisting of 
nuetrons (n), protons (p) and electrons (e). Through the direct (n —> p + e + z7 e ) and/or 
standard (n + n — > n + p + e + v e ) URCA processes the (3 equilibrium 

n <-> p + e~ + v (7) 

can be established. The equilibrium condition requires that the respective chemical poten- 
tials to satisfy 

/in = /ip + /ie + /ip, (8) 

and the charge neutrality requires that 

P P = Pe- (9) 

Since neutrinos do not accumulate in neutron stars a few seconds after their birth fip = 0. 
Within the Fermi gas model one then finds from Eqs. [| and |9] the proton fraction xp = 
PpI (Pn + Pp) hi neutron stars is determined by |39| 

fic^pxp) 1 ^ = AE sym (p)(l - 2x p ). (10) 
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The equilibrium proton fraction is therefore entirely determined by the E sym (p). The values 
of Xp corresponding to the two chosen forms of symmetry energy are shown in the lower 
window of Fig. 1. With the E b ym (p), the xp is zero for p/p > 3, indicating that a pure 
neutron matter could become most stable, leading to the isospin separation instability in HD 
neutron-rich matter. The value of Xp influences many other properties of neutrons stars, 
such as the cooling rate of protoneutron stars [p9|1 . The fast cooling mechanism through 
the direct URCA process can happen if xp is larger than about 1/9. It is seen that with 
the Eg ym (p), the neutron star becomes so proton-rich that the fast cooling can happen at 
densities higher than about 2.3po- On the contrary, it is impossible for this process to happen 
with the E b sym (p). 



C. Equation of state of neutron-rich matter 

With the two forms of the symmetry energy E^ ym (p) and E b ym (p), the corresponding EOS 
for neutron- rich matter is rather different. To illustrate this point, the simplest, momentum- 
independent parameterization 

e( P ,0) = \ + -^u°+ 3 -e%u^ (11) 
I 1 + 0" 5 

is used as the isoscalar part of the EOS, where u = p/po is the reduced density and e° F = 36 
MeV is the Fermi energy. The parameters a = —358.1 MeV, b = 304.8 MeV and o = 7/6 
are determined by saturation properties and a compressibility Koo = 201 MeV of isospin 
symmetric nuclear matter. Our conclusions in this work are qualitatively independent of 
the particular form of e(p, 0). The use of a stiffer EOS, such as = 380 MeV, will 
reduce the compression by about 15% in the energy range studied in this work. This will 
then simply shift the observed isospin effects to slightly higher beam energies. Shown in 
Fig. 2 are the EOS of isospin asymmetric nuclear matter with the E^ ym (p) (upper window) 
and E b sym (p) (lower window), respectively. With the linearly increasing E^ ym (p), the EOS 
becomes stiffer with the increasing 5. The isospin symmetric nuclear matter remains to be 
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the ground state at all densities. This is in stark contrast to the situation using the E b sym (p). 
The EOS obtained with the latter is softened instead of being stiffened by the increasing 
isospin asymmetry 8 at densities higher than 3po- At these high densities, in agreement with 
the information extracted from Fig. 1, the pure neutron matter becomes the most stable 
state whereas the isospin symmetric nuclear matter becomes the most unstable state of HD 
nuclear matter. 

D. Nuclear symmetry potentials and their softest lines 

In accordance with the EOS discussed above, the isoscalar potential is 

v (p) = au + bu a (12) 



and the symmetry potentials can be obtained from 



asy " d Pq 



d Pq " ' 

where q=neutrons and q=protons, respectively. The symmetry potential energy density 

W asy is 

W asy = V 2P 5 2 = [E sym (p) - ci« 2/3 ] P 8\ (14) 

where a = 3/5e^(2 2 / 3 - 1) « 12.7 MeV. With the E a sym (p) one obtains 

< sy = ±2(E sym (p ) - Cl )u5, (15) 

where "+" and "-" signs are for neutrons and protons, respectively. With the E b sym (p), the 
symmetry potential is 

?4 y = ±2(c 2 u c u - c 2 u 2 - Cl u 2/3 )5 + (1 /3cim 2/3 - c 2 u 2 )5 2 , (16) 

where c 2 = E sym (p ) / \u c — 1) w 15 MeV. At low densities v a asy and v h asy are very close since 

limu* w ±2 \c 2 (u c - u) - cim~ 1/3 1 u5 (17) 

U <1 asy l \ J j v / 

w ±2 [c 2 (u c — 1) — ci] u5 

= ±2(E sym (p ) - Cl )uS = v a asy . 
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This is so because of our choice of the two forms of the symmetry energy with the desire to 
study its HD behaviour. This feature is further illustrated in Fig. 3 where the above two 
symmetry potentials are shown as a function of density for 5 = 0.2 and u c = 3. As one 
expects the symmetry potentials are rather similar at p < po but very different at higher 
densities for both neutrons and protons. The different HD behaviour is expected to cause 
neutrons and protons to behave differently during the dynamical evolution of heavy-ion 
collisions. 

Our main purpose in the following will be searching for signs of the different symmetry 
potentials in experimental observables of high energy heavy-ion collisions by comparing 
results for neutrons and protons and using different symmetry potentials. Of course, what 
is more important for the reaction dynamics is the force, or the slope of potentials. With 
the E^ ym (p) the repulsive (attractive) symmetry potential for neutrons (protons) increases 
linearly. With the E b ym (p), however, the symmetry potential for neutrons (protons) has a 
broad maximum (minimum) around p = 1.2p (p = 1.4p ). Thus for both neutrons and 
protons the forces due to the symmetry potentials are zero around the extreme of their 
respective symmetry potentials. This property then leads to a softest point in the resultant 
EOS at a given 5 for isospin asymmetric nuclear matter. Generally, for a varying S there is 
a softest line for both neutrons and protons. These lines are determined by the condition 

(dv b asy /dp) s = 0. (18) 

With the E b sym (p) these lines are along 

, 6ci + 18c 2 m 4 / 3 - 9c 2M( y / 3 

^ftest ~ ± ^ _ Qc2m4/ 3 , (19) 

where the ± sign is for neutrons/protons. These softest lines in neutron-rich matter are 
shown in Fig. 4. In heavy-ion collisions with a proper combination of the beam energy 
and isospin asymmetry, the system can cross these softest lines. The crossing points may 
manifest themselves as local minima in the excitation functions of collective flow for both 
neutrons and protons. This expectation is analogous to that for the local minimum in the 
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excitation function of transverse flow due to the Quark- Gluon-Plasma induced softening of 
the EOS in ultra-relativistic heavy-ion collisions ]5l| . The Coulomb potential may shift the 
relative values and slopes of proton potentials depending on the dynamical evolution of the 
charge distributions during the reaction. It is just difficult to predict quantitatively where 
the minima may appear based on the discussions of the symmetry potentials. Fortunately, 
transport models provide a reliable tool to make more quantitative predictions possible. 

III. ISOSPIN ASYMMETRY OF DENSE NUCLEAR MATTER FORMED IN 
HIGH ENERGY HEAVY-ION COLLISIONS 

High energy heavy-ion collisions provide the only terrestrial situation where the HD 
neutron-rich matter can be formed. Moreover, fast radioactive heavy-ion beams to be avail- 
able at the planned Rare Isotope Accelerator (RIA) in the United States can make the HD 
nuclear matter even more neutron-rich [[51]. The HD behaviour of E sym (p) affects properties 
of the HD nuclear matter formed in high energy heavy-ion collisions. Moreover, interesting 
precursors of the ISI might be observed in collisions with energetic neutron-rich nuclei. 



A. Summary of an isospin-dependent hadronic transport model 

We investigate the effects and phenomena mentioned above within an isospin-dependent 
hadronic transport model for an interacting system of nucleons, Delta resonances and pions 
31| , [52|j53| . Evolutions of the phase-space distribution functions of nucleons, Delta resonances 



and pions with their explicit isospin degrees of freedom are solved numerically by using of the 
test-particle approach |54|,[55|]. Isospin-dependent total and differential cross sections among 
all particles are taken either from the elementary particle scattering data or obtained by 
using the detailed balance [§6|,0 . Explicitly isospin-dependent Pauli blockings for fermions 



are also employed. For a review of the model, we refer the reader to refs. j|,|58| . 

At beam energies above the pion production threshold, baryon resonances play an import 
role in the reaction dynamics. In this study we limit ourselves to beam energies less than 2 
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GeV/nucleon. It is worth mentioning that for medium sized rare isotopes, such as 32 Na, high 
energy beams up to about 950 MeV/nucleon are currently available at GSI and for heavy 
rare isotopes, such as 132 Sn, beams up to 400 MeV/nucleon will be available at the future 
RIA. In this energy range, the most important baryon resonance is the A(1232). Shown in 
Fig. 5 are two examples of the multiplicities of A(1232) resonances and pions. The maximum 
population of the A(1232) is about 2% and 20% in the 132 Sn + 124 Sn reaction at an impact 
parameter of 1 fm and a beam energy of 400 and 2000 MeV/nucleon, respectively. The mean 
field potentials of pions and baryon resonances in nuclear matter are still largely unknown. 
We thus make here a minimum assumption that the isoscalar part of the A potential is the 
same as that for nucleons. To be consistent with the modeling of the isovector potential for 
nucleons, we assume that the isovector potential for A resonances is an average of that for 
neutrons and protons. The weighting factor depending on the charge state of the resonance is 
the square of the Clebsch-Gordon coefficients for isospin coupling in the processes A <-> irN. 
Thus, we have 

v aS y(A-) = v asy (n), (20) 
2 1 

V aS y(A ) = -V aS y(n) + -V aS y(p), (21) 

1 2 

v asy (A + ) = -v asy (n) + -v asy (p), (22) 
v asy (A ++ ) =v asy (p). (23) 

Similarly, we define the effective isospin asymmetry 5; ifce for excited baryonic matter as 

r \Pn)like \Pp)like Ma\ 

d like = 7 s —7 s , (24 ) 

[Pnjlike + [Ppjlike 



where 

2 1 

(Pn)hke = Pn+ gPA° + gPA+ + PA~ , (25) 
2 1 

{Pp)like = Pp+ gPA+ + gPA° + PA++- (26) 

It is evident that the 8uk e reduces naturally to 5 as the beam energy becomes smaller than 
the pion production threshold. 
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B. Isospin asymmetry of dense matter formed in high energy heavy-ion collisions 



Shown in Fig. 6 is the evolution of the central baryon density for the reaction of 
132 Sn + 124 Sn at an impact parameter of 1 fm and beam energies of 200, 400, 1000 and 
2000 MeV/nucleon, respectively. In this energy range a compression of about 1.7 to 3.7 
times the normal nuclear matter density is reached. This very reaction at beam energies 
upto about 400 MeV/nucleon will be available at the RIA facility ||51|. One notices that the 



E sym (p) leads to a slightly higher compression of about 10-15%. This is due to the relatively 
softened nuclear EOS with the E b sym than that with the E^ ym . 

Is the compressed hadronic matter neutron-rich or -poor compared to the initial reaction 
system? The answer to this question has important implications to several critical questions 
in nuclear astrophysics. To answer this question we show in Fig. 7 the correlation between the 
baryon density and the isospin asymmetry 5uk e over the entire reaction volume at the time 
of about the maximum compression in the central 132 Sn + 124 Sn reactions with E/A=400 
(middle window) and 2000 (bottom window) MeV/nucleon, respectively. The whole reaction 
volume is divided into cubic cells of 1 fm 3 in size. Each point in the scatter plot represents 
one such cell. For a comparison, the correlation in the initial state of the reaction is shown 
in the upper window. The initial state for our transport model calculations is generated by 
using neutron and proton density profiles predicted by the Relativistic Mean Field model 



14],|59| . The average 5 of the initial reaction system is 0.22, the low density tail extending to 



the far more neutron-richer side is due to the neutron skin in both the target and projectile 
nuclei. Thus a significant number of cells with densities close to p/po = 1 have isospin 
asymmetries less than 0.22. How this initial p — 5 correlation evolves in a heavy-ion collision 
depends sensitively on the HD behaviour of nuclear symmetry energy. With the E^ ym {p) the 
trend of continuous neutron distillation from higher density regions to lower ones persists at 
all energies. The high density region around p = 2po is about twice more isospin symmetric 
than that with the E b sym (p) at 400 MeV/nucleon. Most interestingly, with the E b ym (p) the 
onset of ISI is clearly seen as indicated by the right turn of 5uk e to higher values at p > 3p 
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in the reaction at Eb eam /A = 2 GeV/nucleon. 

The reason behind the observed neutron-distillation from the HD baryonic matter can 
be easily understood from the density dependence of symmetry energy shown in the upper 
window of Fig. 1. For two parts (1 and 2) of isospin asymmetric nuclear matter to be in 
chemical equilibrium the condition 

E sym (pi)8i = E sym (p 2 )5 2 (27) 

must be satisfied p0|-|63[. Thus, it is energetically favorable in a dynamical process to have 
a migration of nucleons with its direction determined by Eq. |27| according to the density 
dependence of the E sym (p). With the E^ ym (p) which increases linearly with p, there is a 
continuous migration of neutrons (protons) from higher (lower) to lower (higher) densities. 
Whereas with the E b sym (p) which peaks at u = l/2u c , it is most energetically favorable to 
deplete (concentrate) all neutrons (protons) from (to) the peak. Moreover, at u > u c the 
ISI sets in. It leads to an unlimited trapping of neutrons above the critical density u c = 3. 

To be more clear about the roles of the density dependent symmetry energy we investigate 
in Fig. 8 the average (over all cells of the same density) 5uk e as a function of density 
for the central 132 Sn + 124 Sn reactions at 400 and 2000 MeV/nucleon. The overall rise 
of S at low densities is mainly due to the neutron skins of the colliding nuclei and the 
distillated neutrons. Effects due to the different symmetry energies are more clearly revealed 
especially at high densities. For a comparison and study implications of our results to nuclear 
astrophysics, the p — 5 correlation in neutron stars at (3 equilibrium is shown in the insert 
of Fig. 8. The isospin asymmetry at f3 equilibrium 5/3 = 1 — 2xp is completely determined 
by the E sym (p). With the E b ym (p), the 5p is 1 for p/p Q > 3, indicating that the neutron 
star has become a pure neutron matter at these high densities. On the contrary, with the 
E® ym (p), the neutron star becomes more proton-rich as the density increases. An astonishing 
similarity is seen in the resultant 5 — p correlations for the neutron star and the heavy-ion 
collision. In both cases, the symmetry energy E b sym (p) makes the HD nuclear matter more 
neutron-rich than the E® ym (p) and the effect grows with the increasing density. Of course, 
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this is no surprise since the same underlying nuclear EOS is at work in both cases. The 
decreasing E h sym (p) above \u c = 1.5p makes it more energetically favorable to have the 
denser region more neutron-rich. 

The above analyses of the main features are at the instants of about the maximum 
compressions. How do they evolve in time? To answer this question the average n/p ratio of 
the HD region with p / p > 1 is shown as a function of time and beam energy in Fig. 9. The 
effect on (n/p) p / Po >i due to the different E sym (p) is seen to grow with the reaction time until 
the expansion has lead the system to densities below p , especially at higher beam energies. 
Although the compression starts at about the same time, the expansion starts on a faster 
time scale at higher beam energies as one expects. Again, it is seen that whether the HD 
region is neutron-rich or -poor depends critically on the HD behaviour of nuclear symmetry 
energy. 

C. Pion production effects on n/p ratio of dense matter 

During heavy-ion collisions there is a net conversion of neutrons to protons due to the 
more abundant production of 7r~'s than 7r + 's. It is therefore also interesting to explore 
effects of pion production on the n/p ratio of dense matter. This study is also necessary to 
be certain about the effects of the HD nuclear symmetry energy. We thus compare in Fig. 
10 the HD n/p ratio with and without including the pion production channel in the central 
132 Sn + 124 Sn reaction at a beam energy of 2 GeV/nucleon. The beam energy is chosen to 
illustrate the maximum effects of pion production since the latter decreases at lower beam 
energies. We artificially turned off the pion production channel, i.e., N + N — > N + A, and 
attribute its cross section to elastic nucleon-nucleon scatterings such that the total reaction 
cross section is kept a constant. The central densities are shown in the upper rows of Fig. 
10. It is seen that the overall compression is slightly reduced without including the pion 
production. This is because the latter has the role of enhancing the stopping power. It is 
interesting to note that the effect (difference between the curves a and b) on compression due 
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to the different symmetry energy is about the same with or without the pion production. The 
lower two rows demonstrate effects of pion production on the HD n/p ratio. By comparing 
the left and right windows of the lower two rows, it is seen that the pion production reduces 
the HD n/p ratios by about 15%, while the sensitivity to the symmetry energy remains 
almost the same. This finding also indicates that the pion production might be used as a 
perturbative probe of the n/p ratio of HD nuclear matter. 

IV. PROBING THE HD BEHAVIOUR OF NUCLEAR SYMMETRY ENERGY 
WITH HIGH ENERGY HEAVY-ION COLLISIONS 

How to probe experimentally the HD behaviour of E sym (p) in high energy heavy- ion 
collisions? In this section we study four complementary approaches. These are the tt~/tt + 
ratio, nuclear transverse collective flow and its excitation function as well as the neutron- 
proton differential flow. 



A. Pion probe 

As we have shown in the previous section, that the pion production only affects pertur- 
batively the n/p ratio of HD nuclear matter. On the other hand, at beam energies below 
about 2 GeV/nucleon, pions are mostly produced through the decay of A(1232) resonances. 
The primordial tt~/tt + ratio is approximately quadratic in n/p according to the branching 
ratios of single pion production via A resonances in nucleon-nucleon collisions 

--/-+ = - (n/pf. (28) 

op z + np 

Thus the tt~/it + ratio is expected to be a sensitive probe of the n/p ratio of HD nuclear 
matter. Pion reabsorptions and rescatterings (it + N — > A and N + A — ► N + N) are 
expected to complicate the above relationship. Nevertheless, very high sensitivity to the 
n/p ratio is retained in the final 7i~ /tt + ratio as indicated in the experimental data of high 



energy heavy-ion collisions [64|. Shown in Fig. 11 are the (tt /n + )m i . e ratio 
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as a function of time for the central 132 Sn+ 124 Sn reaction at beam energies from 200 to 2000 
MeV/nucleon. This ratio naturally becomes the final tt~ /tt + ratio at the freeze-out when the 
reaction time t is much longer than the lifetime of the delta resonance ta- The (tt~ /7i + )uk e 
ratio is rather high in the early stage of the reaction because of the large numbers of neutron- 
neutron scatterings near the surfaces where the neutron skins of the colliding nuclei overlap. 
By comparing Fig. 9 and Fig. 11, it is seen that a variation of about 30% in the (n/p) P /p >i 
ratio due to the different E sym (p) results in about 15% change in the final tt^/ty + ratio. It 
has thus an appreciable response factor of about 0.5 to the variation of HD n/p ratio and is 
approximately beam energy independent. Therefore, one can conclude that the (n^ /n + )hke 
ratio is a direct probe of the HD n/p ratio, and thus an indirect probe of the HD behaviour 
of nuclear symmetry energy. 



B. Transverse collective flow and its excitation function 

Next we investigate the transverse collective flow as a probe of the HD behaviour of 
nuclear symmetry energy. We perform the standard analysis of the average transverse mo- 
mentum in the reaction plane |^5| 



1 N(y) 

<Px /N>(y) = — Y { PUy)i (30) 

where N(y) is the number of particles at rapidity y and pi X is i th particle's transverse 
momentum in the reaction plane. Shown in Fig. 12 are two typical examples of the transverse 
collective flow analysis for neutrons and protons for the mid-central 132 Sn + 124 Sn reactions 
at 400 MeV/nucleon. Appreciable effects of the nuclear symmetry energy is clearly seen, 
especially for protons. 

We found that the symmetry energy effects on the transverse flow vary with beam energy 
and are different for neutrons and protons. To better characterize the effects we study in 
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Fig. 13 the excitation function of flow parameter. The latter is defined as the slope of the 
transverse momentum distribution at mid-rapidity, i.e. 



We concentrate on analyzing differences between the collective flows of neutrons and protons 
as well as those caused by using the different symmetry energies. These differences can 
be mainly attributed to the different symmetry potentials for neutrons and protons since 
the isoscalar potential has no direct contribution to these differences to the first order of 
approximation. On top of the generally growing flow parameter there are distinct structures 
closely associated with the symmetry energy used. In particular, a local minimum appears in 
the excitation functions for both neutrons and protons around Eb eam /A = 500 MeV with the 
E b sym . The dip is more obvious for protons than neutrons. These minima are direct signals 
of the softening of the underlying EOS as we expected from studying the softest lines of 
the symmetry potentials. With the E a sym there is also a local minimum in the excitation 
function for protons around E beam /A = 300 MeV. This minimum is a result of the interplay 
between repulsive hadronic scatterings and the attractive symmetry potential for protons 
(the solid line in the bottom window of Fig. 3). To be more quantitative, the average number 
of successful hadronic scatterings per nucleon and the maximum baryonic compression are 
shown as a function of beam energy in Fig. 14. Because the E b sym softens the EOS it thus 
leads to an about 3% higher compression and more scatterings compared to the case with 
the E° 1 . The small increase in compression has very little effect on the scalar potential. 
Thus the observed structures can not be due to the isoscalar potentials. Below the pion 
production threshold in nucleon-nucleon collisions at about 296 MeV, the Pauli blocking is 
still strong and thus the collision numbers increase very little with beam energy. Whereas 
above the pion production threshold, because of the opening of the inelastic channel via 
N + N — > N + A and its resonance nature, the collision numbers grow quickly. With the 
E° , for protons the symmetry potential causes negative flow while the collisions result in 
positive flow. Around 300 MeV/nucleon the collision numbers are still small, effects of the 




(31) 
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symmetry potential are thus relatively important in the interplay which leads to the dip near 
the pion production threshold. At higher beam energies, however, the flow increases because 
of the dominating role of the collisions and their fast growth with beam energy. While for 
neutrons with the E® ym , both the collisions and the repulsive symmetry potential (the solid 
line in the upper window of Fig. 3) enhance the transverse flow. Thus the flow parameter 
for neutrons increases almost monotonously with beam energy. Therefore, the structures 
observed above in the excitation functions of the nuclear transverse flow parameter are 
unique qualitative signatures of the density dependence of symmetry energy. 



C. Neutron-proton differential flow probe 



The neutron-proton differential collective flow is measured by |33 

1 N(y) 

F np(y) = E VxJh (32) 

where N(y) is the total number of free nucleons at the rapidity y, p Xi is the transverse 
momentum of particle i in the reaction plane, and r, is +1 and —1 for neutrons and protons, 
respectively. The free nucleons are identified as those having local nucleon densities less than 
l/8po- The F np (y) combines constructively the in-plane transverse momenta generated by 
the isovector potentials while reducing significantly influences of the isoscalar potentials of 
both neutrons and protons. Thus, it can reveal more directly the HD behaviour of E sym (p) in 
high energy heavy-ion collisions. Two typical results for mid-central 132 Sri + 124 Sn reactions 
at 400 MeV/nucleon and 1000 MeV/nucleon, respectively, are shown in Fig. 15. A clear 
signature of the HD behaviour of E sym (p) appears at both forward and backward rapidities. 
To characterize the effect, the slope dF np / d(y cm / ybeam) at mid-rapidity is shown as a function 
of beam energy in Fig. 16. A striking difference of about a factor of 2 exists in the reactions at 
Ebeam > 200 MeV/nucleon. This large effect can be very easily observed by using available 
detectors at several heavy-ion facilities in the world. However, the interesting structures 
observed in the excitation function of the transverse flow parameter F is largely being 
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smeared out because the different rapidity distributions of neutrons and protons also enters 
in the calculation of the neutron-proton differential flow. Compared with the n~ /n + ratio, 
both the transverse flow and the neutron-proton differential collective flow are more directly 
affected by and are thus also more sensitive probes of the HD behaviour of E sym (p). 



V. SUMMARY 



In summary, the HD behaviour of nuclear symmetry energy has been puzzling physi- 
cists for decades. In this work, high energy heavy-ion collisions are proposed as a novel 
means to solve this longstanding problem. Within an isospin dependent hadronic transport 
model using two representative density dependent symmetry energy functions predicted by 
many body theories, it is shown that the isospin asymmetry of HD nuclear matter formed 
in high energy heavy-ion collisions is uniquely determined by the HD behaviour of E sym (p). 
It is demonstrated that the isospin separation instability indeed can happen in high energy 
heavy-ion collisions. We explored several promising experimental probes for the HD symme- 
try energy. Among them, four quantitatively sensitive observables, the tc~ /n + ratio, nuclear 
transverse collective flow and its excitation function as well as the neutron-proton differ- 
ential collective flow are studied. A precursor of the possible isospin separation instability 
in dense neutron-rich matter is predicted to appear as the local minima in the excitation 
functions of the flow parameter for both neutrons and protons above the pion production 
threshold. This precursor can be used as a unique signature of the isospin dependence of the 
nuclear EOS. This qualitative signature is advantageous over the other quantitative observ- 
ables because the latter are often affected by several ingredients of the reaction dynamics 
and experimental uncertainties. It should be noted that the present study is based on a 
momentum-independent transport model. The momentum-dependent potential is known to 
be important for a reliable description of global momentum distributions. Nevertheless, we 
expect that all qualitative features, such as the existence of the minimum in the excitation 
function of transerse flow, are not affected by the momentum-dependence. However, the 
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quantitative results, such as the exact location of the flow minimum, are affected by the 
momentum-dependence of the nuclear potential. A more refined study using a momentum- 
dependent isoscalar potential is underway and the results will be reported elsewhere. It 
should also be mentioned that the fragmentation mechanism of isospin asymmetric nuclear 
matter is still not clear and remains a hot topic of current studies . The different frag- 
mentation mechanisms might affect the neutron/proton ratio of the observed free nucleons, 
and thus the neutron-proton differential flow. It would thus be useful to measure also ob- 
servables associated with fragments together with the single particle observables examined in 
this study. It still remains a theoretical challenge to develop a fully quantum transport the- 
ory that deals with both particle production and cluster formation in a numerically tractable 
manner. The present work shows the need and promise of working in this direction. Ex- 
perimental measurements of the proposed observables will provide the first terrestrial data 
to constrain stringently the HD behaviour of nuclear symmetry energy. Future compar- 
isons between the experimental data and model calculations will allow us to extract crucial 
information about the EOS of dense neutron-rich matter. 

This work was supported in part by the National Science Foundation Grant No. PHY- 
0088934 and Arkansas Science and Technology Authority Grant No. 00-B-14. 
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FIGURES 




FIG. 1. Upper window: Two representatives of the nuclear symmetry energy as a function of 
density. Lower window: the corresponding proton fractions in neutron stars at (3 equilibrium. 
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FIG. 4. The softest lines for neutrons and protons using the symmetry potential corresponding 
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FIG. 5. Multiplicities of Delta resonances and pions in the central reaction of Sn+ Sn at 
a beam energy of 400 MeV/nucleon (upper window) and 2 GeV/nucleon (lower window). 



30 




t (fm/c) 



FIG. 6. Evolution of the central baryon density in the central reaction of 132 Sn + 2 Sn at a 
beam energy from 200 to 2000 MeV/nucleon. 
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FIG. 7. Scatter plots of isospin asymmetry-density correlation in the initial state (upper win- 
dow) and at the instants of maximum compressions in the reaction of 132 Sn + 124 Sn with the 
nuclear symmetry energy E% (filled circles) and E\ \ (right triangles), respectively. 



32 




FIG. 8. Upper window: the isospin asymmetry-density correlations at t=20 fm/c and 
Ebeam/A = 400 MeV in the central 132 Sn + 124 Sn reaction with the nuclear symmetry energy 
E a sym and E b sym , respectively. Lower window: the same correlation as in the upper window but at 
10 fm/c and E^ eam /A = 2 GeV/nucleon. The corresponding correlation in neutron stars is shown 
in the insert. 
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FIG. 9. The neutron/proton ratio of nuclear matter with density higher than the normal nuclear 



matter density as a function of time with the nuclear symmetry energy E a sym and E b sym , respectively. 
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FIG. 10. Central baryon density (top row) and the neutron/proton ratio of nuclear matter 
with density p > 2po (middle row) and p > po (bottom row) with (right colomn) and without (left 
colomn) the pion production channel in the central 132 Sn + 124 Sn reaction at a beam energy of 2 
GeV/nucleon. 
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FIG. 11. 7r /tt + ratio as a function of time in the central 132 Sn + 124 Sn reaction at a beam 
energy of between 200 and 2000 MeV/nucleon. 
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FIG. 12. Transverse flow analysis for nuetrons (upper window) and protons (lower window) 
in the mid-central 132 Sn + 124 Sn reaction at a beam energy of 400 MeV/nucleon. 
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FIG. 13. The excitation functions of the transverse flow parameter for neutrons (upper) and 
protons (lower) in the mid-central 132 Sn+ 124 Sn reactions. 
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FIG. 14. The average number of successful hadronic scatterings per initial nucleon (upper) and 
the maximum baryon density reached in the mid-central 1? ' 2 Sn + 124 Sn reaction as a function of 
beam energy. 
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FIG. 15. The neutron-proton differential collective flow in the mid-central 132 Sn + 124 Sn re- 
actions at Eb eam /A = 400 MeV (upper window) and 1000 MeV (lower window) with the nuclear 
symmetry energy E a sym and E b sym , respectively. 
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FIG. 16. Excitation function of the slope parameter of the neutron-proton differential flow for 
the mid-central 132 Sn + 124 Sn reaction. The lines are drawn to guide the eye 
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